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Abstract. Wc discuss the relaxation and dephasing rates that result from the control and the measurement 
setup itself in experiments on Josephson persistent-current qubits. For control and measurement of the 
qubit state, the qubit is inductively coupled to electromagnetic circuitry. We show how this system can be 
mapped on the spin-boson model, and how the spectral density of the bosonic bath can be derived from 
the electromagnetic impedance that is coupled to the qubit. Part of the electromagnetic environment is 
a measurement apparatus (DC-SQUID), that is permanently coupled to the single quantum system that 
is studied. Since there is an obvious conflict between long coherence times and an efficient measurement 
scheme, the measurement process is analyzed in detail for different measurement schemes. We show, that 
the coupling of the measurement apparatus to the qubit can be controlled in situ. Parameters that can be 
realized in experiments today are used for a quantitative evaluation, and it is shown that the relaxation and 
dephasing rates that are induced by the measurement setup can be made low enough for a time-resolved 
study of the quantum dynamics of Josephson persistent-current qubits. Our results can be generalized as 
^ ■ engineering rules for the read-out of related qubit systems. 

PACS. PACS-03.67.Lx Quantum computation - PACS-05.40.-a Fluctuation phenomena, random pro- 
cesses, noise, and Brownian motion - PACS-74.50.+r Proximity effects, weak links, tunneling phenomena, 
and Josephson effects - PACS-85.25.Dq Superconducting quantum interference devices (SQUIDs) 

, 1 Introduction it was recognized that the validity of quantum mechanics 

' f° r a macroscopic degree of freedom could be tested in real 

The dynamics of electromagnetic circuits and other macro- experiments Q . 

I ' scopic objects is usually well described by classical laws: T T . , , ,. ... 

*7-j ■ , i , , ,., ... In 1980, Lcggctt pointed out that cryogenic and mi- 

V , quantum coherent phenomena like superposition states . . . , .' , : , . . . , , , , , 

CZ n a u j ■ a rru c j crolabrication technologies had advanced to a level where 

are usually not observed m macroscopic systems, the rounders .„,..,. ° , , , ., , , 

r , i , u • i j-i i j i i-i i macroscopic bchrodmger s cat states could possibly be re- 

Ool the quantum mechanical theory already recognized mat ,. i • ,, , . , , • T i 

. . r , n- . i , , • -I, r j ahzed m small superconducting loops that contain J oscph- 

• • ■ there is m tact a conflict between a straight forward ex- . . . .r-iQ, T , , , , 

l>,. i i • p , , . , . , son tunnel mnctions BU4 . In such systems, the Joseph- 

trapolation ol quantum mechanics to a macroscopic scale, . / . . Lt - t v l . . ' . 

, , r i i i ■ ,i . ,i son phase or equivalently. the persistent supercurrcnt in 

and the laws oi classical physics that govern the macro- . f . > . ,. r ^ n 

ii x ,. ! X ■ -u-i-4- c the loop) is a collective coordinate for the Cooper-pair 

scopic world. In particular, this concerns the possibility ol , • • . , . , 

r-t , ... r ,1 .. j. j. c . condensate, and it is conjugate to a variable which de- 

quantum superpositions ol collective coordinates (i. e. center- .. ' J ° T 

of-mass-like coordinates) of objects that are much bigger sc \ lhes the ch ^ r S e difference across the Josephson tun- 

than the atomic scale. These difficulties were first pre- nel J unctl0n - However, while the analysis of the isolated 
sented by Schrodinger @, and are now known as Schr6dinger'& uantum s ^ tem shows that superpositions of the macro- 

t. i o u "A' j j' • fiv. i • ji i scopic coordinates might very well occur m these loops, 

cat paradox. Schrodinger s discussion ol the cat m the box ... . . , , , , • , , 

i , , , , i /~\ i i it is by no means obvious that such behavior can also be 
was clearly meant as a qeaanken experiment. Only several , , , . „ i • • 
decades later, after the discovery of the Josephson effect, demonstrated experimentally. Such superposition states 
are extremely iragile, reflecting the tendency oi macro- 
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but permanently placed in an experimental setup may 
hinder attempts to study macroscopic quantum coher- 
ence. Nevertheless, interesting results with evidence for 
macroscopic quantum tunneling, energy level quantiza- 
tion and coherent dynamics between quantum levels were 
obtained with systems where the Josephson phase coor- 
dinate is trapped in a metastable well (for an overview 
see |5|pMP ,Fo| ) ■ Also with systems where the energy 
scale for single-charge effects is higher than, or compa- 
rable to the energy scale for the Josephson effect, quan- 
tum coherent dynamics has recently been demonstrated 
|jlf|,^2). In Josephson junction loops, quantum superposi- 
tion states of persistent currents have been demonstrated 
spectroscopically [ p^JT^I . However, time-resolved experi- 
ments that prove quantum-coherent oscillations between 
macroscopically-distinct persistent-currents states in the 
sense of Ref. jl have not been reported yet. 

Whether such experiments can be realized at all has 
been intensively discussed in the literature flUl , without 
consensus being reached. However, a detailed analysis with 
estimates based on measurement techniques that can be 
realized in experiments today, has been discussed very lit- 
tle. The quantum coherent dynamics observed with the 
other Josephson junction systems (such as Cooper pair 
boxes [ll],[l2| or single junctions [ppc||), indicates that it 
might be possible to obtain similar experimental results 
with Josephson persistent-current loops. Efforts in this 
direction were stimulated by the prospect that it might 
be possible to realize a quantum computer with supercon- 
ducting Josephson devices |T|,|l^,^8|||g|l| • An impor- 
tant advantage of a Josephson quantum computer would 
be that, if accurate quantum coherent control of elemen- 
tary units would be possible, it would be a system that 
can be extended to one containing a very large number of 
quantum bits (qubits). The large size of the qubits allows 
for individual (local) control and readout of the qubits and 
qubit-qubit couplings. 

In this article we analyze the feasibility of demonstrat- 
ing quantum coherent dynamics of Josephson persistent 
currents with experimental techniques for manipulating 
and reading qubit states that can be realized in the lab- 
oratories today (i. e. assuming the available techniques 
for device fabrication, cryogenics, microwave applications 
and electronic filtering). Such mesoscopic solid-state ex- 
periments suffer from the fundamental difficulty that one 
cannot avoid that an electronic measuring device is per- 
manently coupled to the single quantum system that is 
studied We will not consider future measurement 

techniques which may couple less directly to the qubit. 
A meter must be present in any useful experiment, and, 
unlike experiments with for instance photons, this means 
that a measuring device must be permanently located very 
close to the solid-state quantum bit (e. g. fabricated on the 
same chip). With such a setup, there is obviously a conflict 
between an efficient measurement scheme with a strong 
measurement, and long decoherence times in the quan- 
tum system that is studied. For successful experiments in 
this direction, a detailed understanding of the measure- 
ment scheme is therefore needed such that the decoher- 



ence that is induced by the setup itself can be reduced to 
an acceptable level. 

Obviously, there exist many other sources of decoher- 
ence for Josephson qubits that one should worry about as 
well. The critical current of the junctions may show tele- 
graph noise p3| , which would give rise to decoherence sim- 
ilar to what is described in Ref. [Q. Moreover, it has been 
stressed that a very high number of spin degrees of free- 
dom is usually present in the solid state environment that 
may decohere Josephson qubits (see the work by Prokof 'ev 
and Stamp [^5|,^6) on the spin-bath, and Ref. |27j] for esti- 
mates for persistent-current qubits). Another example is 
decoherence from quasiparticles that effectively shunt the 
junction ^8|. These effects themselves are very interesting 
for further study. However, a study of for example the de- 
phasing due to spin impurities remains impossible as long 
a reliable and well-understood measurement scheme for 
the loop's quantum dynamics is not available. Therefore, 
we will concentrate here on dephasing and mixing due to 
the experimental wiring and the measurement scheme it- 
self. 

Our analysis mainly focuses on experiments with the 
three-junction persistent-current qubit proposed by Mooij 
et at, [p^^oUl^] , in a setup where they are measured by 
undcrdamped DC-SQUID magnetometers (in this article 
we will reserve the word SQUID for the measuring DC- 
SQUID (Fig. |l|a), and not use it for the three-junction 
qubit (Fig. |], center)). The decohering influence of the 
inductively coupled DC-SQUID is analyzed as well as de- 
coherence that results from inductive coupling to on-chip 
control lines for applying microwave signals and local mag- 
netic fields. Model descriptions of the experimental setup 
will be mapped on the spin-boson model, such that we 
can use expressions for the relaxation and dephasing rates 
from the spin-boson literature. The typical experimental 
situations will be described quite extensively to justify the 
models and the approximations used. The results will be 
worked out quantitatively, and we will evaluate whether 
we can realize mixing and dephasing rates that are com- 
patible with measurement schemes based on DC-SQUIDs. 
The design criteria developed in this work are more gen- 
eral and should also be of interest for experiments on loops 
with a single Josephson junction fll3|] , and quantum cir- 
cuits where the charge degree of freedom is measured, as 
Josephson charge quantum bits |l8[|2l|| and quantum dots 
p9[ . In a more general context the value of this work is 
that it presents in detail an example of a measurement 
process on a single quantum system in which the deco- 
herence enhances with increasing measurement strength. 
The issues discussed here are an example of experimen- 
tal difficulties that will unavoidably play a role in many 
realizations of quantum computers. 



1.1 Outline 

In section 2 we will summarize a theoretical description 
of the Josephson persistent-current qubit, and the spin- 
boson theory that will be applied in our analysis. Section 
3 presents a description of the measurement process with 
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the DC-SQUIDs, and a typical scheme for coupling the 
qubit to the on-chip control lines. In section 4 we work out 
the qubit 's relaxation and dephasing rate that result from 
the coupling to a switching DC-SQUID. This is worked out 
quite extensively and the definitions presented in this sec- 
tion are also used in section 5. Two measurement scenarios 
with different types of electromagnetic shunt circuits for 
the DC-SQUID will be compared. A short analysis of the 
decoherence due to the coupling to on-chip control lines 
is presented in section 5. Section 6 presents a few control 
techniques that can improve decoherence rates. 



2 Qubit Hamiltonian and theory for 
relaxation and dephasing 

This work aims at calculating relaxation (mixing) rates 
and dephasing (decoherence) rates for a Josephson persistent 
current qubit which result from its inductive coupling to 
the measurement setup. The measurement setup is formed 
by a DC-SQUID and control lines, which are attached 
to leads and coupled to filters and electronics (Fig. |f|). 
This setup will be modeled as a macroscopic quantum 
two-level system (central spin) that is coupled to a linear 
electromagnetic impedance Z t (uj), where lo the angular 
frequency. The impedance Z t (uj) forms an oscillator bath 
and can be described by a set of LC oscillators. This al- 
lows for mapping the problem on the spin-boson model: 
a central spin-i system that is coupled to a bosonic bath 
|[(||l| . The parameters of the bath will be derived from 
the Johnson-Nyquist noise from Z t {uj). In this section we 
will first introduce the qubit Hamiltonian and physical 
properties of the qubit, and then summarize the spin- 
boson expressions for relaxation and dephasing. 



2.1 Qubit properties and Hamiltonian 



A as the tunnel splitting at <I> q — j^q, such that A = 
2W with W the tunnel coupling between the persistent- 
current states. This system has two energy eigen values 
±-|v 'A 2 + e 2 , such that the level separation v gives 



v= y/A 2 



(2) 



The three- Josephson junction qubit fl 9| , |20| , |l 4| is a low- 
inductance superconducting loop which contains three Joseph- 
son tunnel junctions (Fig. [I]). By applying an external 



In general A is a function of s. However, it varies on 
the scale of the single junction plasma frequency, which 
is much above the typical energy range at which the qubit 
is operated, such that we can assume A to be constant for 
the purpose of this paper. 

In the experiments <P q can be controlled by applying a 
magnetic field with a large superconducting coil at a large 
distance from the qubit, but for local control one can apply 
currents to superconducting control lines, fabricated on- 
chip in the direct vicinity of the qubit. The qubit's quan- 
tum dynamics will be controlled with resonant microwave 
pulses (i. e. by Rabi oscillations). The proposed opera- 
tion point is at e ~ 5A, which was analyzed to be a good 
trade-off between a system with significant tunneling, and 
a system with o^-like eigen states that can be used for 
qubit-qubit couplings and measuring qubit states ||l9| , po| . 
For optimal microwave control the qubit will be placed in 
a small off-resonant cavity, and the microwave signals will 
be applied through on-chip superconducting control lines 
(i. e. the magnetic component of the fields from microwave 
currents will be used). The qubit has a magnetic dipole 
moment as a result of the clockwise or counter-clockwise 
persistent current The corresponding flux in the loop is 
much smaller than the applied flux <P q , but large enough 
to be detected with a SQUID. This will be used for mea- 
suring the qubit states. For our two-level system Eq. [l], 
this means that both manipulation and readout couple to 
a z . Consequently, the noise produced by the necessary cir- 
cuitry will couple in as flux noise and hence couple to a z , 
giving e a small, stochastically time-dependent part 5e(t). 
Our system also has electric dipole moments, represented 
by o x . These couple much less to the circuitry and will 
hence not be discussed here. 



flux $ q a persistent supercurrent can be induced in the 
loop. For values where <P q is close to a half-integer num- 
ber of superconducting flux quantums $o, two states with 
persistent currents of opposite sign are nearly degenerate 
but separated by an energy barrier. We will assume here 
that the system is operated near <P q — |^o- Classically, 
the persistent currents have here a magnitude I p . Tunnel- 
ing through the barrier causes a weak coupling between 
the two states, and at low energies the loop can be de- 
scribed by a Hamiltonian in the form of a two-level system 

PPPL 

-77°x, (1) 



2.2 Spin-boson theory for relaxation and dephasing 

For defining the relaxation and dephasing rates, the state 
of the qubit is described with a reduced density matrix p, 
in the basis which is spanned by the eigen vectors of <r z in 
(|l|), i.e. by the semiclassical states with well-defined left 
(L) or right (R) circulating current 



P : 



Pl,l Pr,l 

PLM PR,R 



(3) 



H q = \o z 



where a z and a x are Pauli spin operators. The two eigen 
vectors of a z correspond to states that have a left or a 
right circulating current and will be denoted as \L) and 
\R). The energy bias e = 2I p (<P q — g$o) is controlled by 
the externally applied field S q . We follow p3] and define 



We will concentrate our discussion on the undriven case. 
The qubit dynamics consists of quantum-coherent oscilla- 
tions, which decay on a time-scale = , the dephas- 
ing time. This dephasing is superimposed on an energy 
relaxation mechanism on a larger timescale r r = r~ , 
the relaxation time. This combined decoherence process 
brings the system into an incoherent thermal mixture of its 
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energy eigen states. Expressed in the basis of these eigen 
states, the off-diagonal terms (coherences) of the density 
matrix p go to zero on the time scale of , whereas the di- 
agonal terms (populations) decay in r r to the Boltzmann 
factors. For estimating r r and we will work from the 
systematic weak-dam ping approximation (SWDA) devel- 
oped by Grifoni et al. (32] , which covers recent theoretical 
progress for the spin-boson theory. Grifoni et al. calculated 
expressions for r r and r<f, for a spin-boson system in which 
the coupling to the environment is dominated by bilinear 
coupling terms between a z and the bath coordinates. This 
is a good approximation for a quantum two-level system 
that is only weakly damped by the environment. 

In our case the bath is formed by the impedance Z t (oj), 
and can be described by a set of LC oscillators with flux 
coordinates conjugate charge coordinates Qi, and Hamil- 
tonian 

Hbath = ^ ($1/ 2L * + Q?/ 2 Ci) ■ ( 4 ) 

The flux produced by the qubit will shift the flux <Pi in 
each LC oscillator. The coupling Hamiltonian is 



which, if J(u>) is a sufficiently smooth function of lo can 



usually be taken as a 



1 dJ(w) 



at lo sa 0. These results 



H, 



(j z \ ^ » 

q-bath = — / ,. Ci^i 



2 



(5) 



where c, is the coupling strength to the i-th oscillator. 
In this model the influence of the oscillator bath on the 
qubit can be captured in the environmental spectral den- 
sity function 



2tt du 

can be intuitively interpreted: The system can relax by 
dissipating all its energy v into an environmental boson. 
Due to the weakness of the coupling, there are no multi- 
boson processes. The relaxation also dephases the state. 
Moreover, dephasing can occur due to the coupling to low- 
frequency modes which do not change the energy of the 
system. These expressions for relaxation and dephasing 
have also been found by studying the Hamiltonian of our 
qubit coupled to a damped oscillator, usin g a Markovian 
master equation approach by Tian et al. fl34f (based on 
work by Garg et al. J3f|). 

The expressions (0) and (||) have prefactors (77) 2 and 

(j) that depend on the tunnel splitting A and the energy 
bias e. These factors correspond to the angles between 
noise and eigen states usually introduced in NMR |7J| and 
account for the effect that the qubit's magnetic dipole ra- 
diation is strongest where the flux in the qubit <P q = ^<Po 

(i. e. (^7) maximal), and that the level separation v is in- 
sensitive to flux noise at this point (i. e. ^7 = (^) w 0). 
One should know and control J(lo) at the frequency v/K 
for controlling the relaxation, and at low frequencies for 
controlling the dephasing. In this article we will calculate 
the noise properties of a few typical experimental environ- 
ments, and calculate how the noise couples to the qubit. 
This can be used to define J(u>) for our specific environ- 
ments. 



(6) 



where LOi the resonance frequency of the i-th oscillator. 
The dense spectrum of the degrees of freedom in the elec- 
tromagnetic environment allows for treating J(lo) as a con- 
tinuous function. 

From now on, we focus on the low-damping limit, J(u>) -C 
lo. Thus, the energy-eigenstates of the qubit Hamiltonian, 
Eq. [l], are the appropriate starting point of our discussion. 
In this case, the relaxation rate r r (and relaxation time 
r r ) are determined by the environmental spectral density 
J{lo) at the frequency of the level separation v of the qubit 



3 Measurement setup 

This section describes a typical experimental setup for 
measurements on Josephson persistent-current qubits. Only 
the parts that are most strongly coupled to the qubit 
will be worked out (Fig. |l|). The first part describes a 
DC-SQUID magnetometer that is used by measuring its 
switching current, the second part addresses the use of on- 
chip superconducting lines for applying magnetic fields to 
the qubit. 



r r = t- 1 = 



2 V v ) 1 1 ' \2k B T 



(7) 



where T is the temperature of the bath. The dephasing 
rate (and dephasing time t^,) is 



t et, = T, = — 



-] a2u- 

v i 



k R T 



(8) 



These expressions have been derived in the context of 
NMR p3| using a Markov approximation and recently 
been connrmed by a full path-inegral analysis p^ j. 

The second term only contributes for an environment 
which is Ohmic at low frequencies (i. e. for J(uj) cx lo). 
Here a is a dimensionless dissipation parameter. It is de- 
termined by the slope of J(lo) at low frequencies 



= lim 



J(w) 



+0 2-KLO 



(9) 



3.1 Switching DC-SQUID 

SQUIDs are the most sensitive magnetometers, and they 
can be operated at very low power consumption |36| . We 
will consider here the use of a DC-SQUID with a hys- 
teretic current-voltage characteristic (IV), and unshunted 
junctions that are extremely underdamped. It is used by 
ramping a current through it and recording the switch- 
ing current: the bias current at which it switches from 
the supercurrent branch to a nonzero voltage in its IV 
(Fig. ^). The switching current is a measure for the mag- 
netic flux in the loop of the SQUID. An important ad- 
vantage of this scheme is that the SQUID produces before 
readout very little noise. As long as the SQUID is on the 
supercurrent branch, it does not produce any shot noise or 
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Josephson oscillations. If the external noise and interfer- 
ence can be suppressed by filtering, there is only Johnson- 
Nyquist noise from the low-temperature leads and filtering 
that the SQUID is connected to. At low frequencies this 
residual noise has little power since the device is supercon- 
ducting. Moreover, we will show in section 4 that at low 
bias currents the effective coupling between this meter and 
the quantum system is very weak. In comparison, damped 
non-hysteretic SQUIDs have the problem that the shunt 
resistors at the junctions also provide a damping mecha- 
nism for the qubit. In a hysteretic SQUID there is more 
freedom to engineer the effective impedance seen by the 
qubit, and it also has the advantage that the voltage jump 
at the switching current is much larger |37]|. Recently, a 
similar scheme with a superconducting single-charge de- 
vice, that can be operated as a switching electrometer has 
been reported |58],[l2) • Voltage biased single-electron tran- 
sistors fo r q uantum measurements have been analyzed in 
Refs. [|9|gol0,||. 

For qualitative insight in the measurement process we 
will present here a simplified description of the SQUIDs 
noise and dynamics (valid for a DC-SQUID with symmet- 
ric junctions and a loop with negligible self inductance). 
In section 4 it will be worked out in more detail. The su- 
percurrent through the SQUID with a flux <P in its loop 
is 



cos/ sinipe 



(10) 



where / = ir<I>/<Po, I co the critical current of the junc- 
tions, and (p ex t a Josephson phase coordinate. I sq will 
be distinguished from the applied bias current Iuas, as 
part of the bias current may go into circuitry shunting the 
SQUID. Insight in the SQUID's response to a bias current 
is achieved by recognizing that ( |io| ) gives steady state so- 
lutions (dU I d(p ex t = 0) for a particle with coordinate ip ex t, 
trapped on a tilted washboard potential (Fig. ^) 



U 



2c 



(2I co COS / COS (p ex t + I s q<Pext) ■ (H) 



In this picture, the average slope of the potential is propor- 
tional to the bias current, and the supercurrent branch of 
the SQUID's IV corresponds to the particle being trapped 
in a well. The Josephson voltage across the SQUID V = 
^k~it L * s nonzero for the particle in a running mode. In 
absence of noise and fluctuations, the SQUID will switch 
to the running mode at the critical current Ic 



Ir 



21 co | cos /| 



(12) 



A DC-SQUID can thus be regarded as a single Joseph- 
son junction with a flux-tunable critical current. In prac- 
tice, noise and fluctuations of <p ex t will cause the SQUID 
to switch before the bias current reaches Ic- This cur- 
rent level will be denoted as the switching current Isw 
to distinguish it from Ic- It is a stochastic variable, but 
averaging over repeated recordings of Isw allows for de- 
termining / with great accuracy. This naive description 
can be used to illustrate three important properties of the 
measurement process with the SQUID. 



In the experiment, the electronics for recording the 
SQUID's IV obtains information about / when the SQUID 
switches. However, rewriting ( |l0| ) as 



<Pext 



sm 



2I co cos / 



(13) 



shows that the SQUID's coordinate (p ex t is already corre- 
lated (i. e. entangled) with the flux / at current values be- 
low Isw- Small voltage fluctuations that result from small 
plasma oscillations and translations of <p e xt will cause dis- 
sipation in the electromagnetic environment of the SQUID, 
which damps the dynamics of ip ex t- This means, that in 
a quantum mechanical sense, the position of (p ex t, and 
thereby /, is measured by the degrees of freedom that 
form the electromagnetic impedance that is shunting the 
SQUID (i. e. the leads and filtering between the SQUID 
and the readout electronics), and that the measurement 
may in fact take place before it is recorded by a switching 
event. 

Secondly, ( |l5| ) shows that the SQUID's coordinate <p ex t 
is independent of the flux in the loop (dip ex t/df = 0) 
for I sq = 0. Therefore, in absence of fluctuations of <p ex t 
and current noise, the meter is at zero current effectively 
"off" . In practice this can not be perfectly realized, but it 
illustrates that the decoherence from the SQUID may be 
reduced by a large extent at low bias currents. 

Thirdly, for bias currents well below Ic, the coordinate 
ip e xt is trapped in a potential that is for small oscillations 
close to harmonic. The SQUID can in this case be regarded 
as an inductance 



Lj = 



1 



2e 



4&COB 2 / - P sa 



(14) 



(see also ( |26| ) below). The noise from the SQUID can 
here be described by the Johnson- Nyquist noise from the 
SQUID's Josephson inductance jl4[ ) in parallel with the 
SQUID's environmental impedance (Fig. ||a, b). For high 
bias currents very close to Ic, the spectrum will have more 
power and calculating the noise properties will be more 
complicated. Here non-harmonic terms in the trapping 
potential become important, and there maybe additional 
noise from a diffusive motion of ip ex t to neighboring wells 
(Fig. |c). For hysteretic SQUIDs this regime with diffu- 
sive motion of ip ex t and switching currents very close to 
Ic will only occur in SQUIDs with a very specific electro- 
magnetic shunt . In many realizations of hysteretic 
DC-SQUIDs ip e xt will escape to a running mode without 
retrapping in lower wells (Fig. |^d), and Isw can be much 
lower than Ic- In this case the approximation using ( |l4| ) 
should be valid for description of the noise before a switch- 
ing event. 

The statistics of Isw readouts depend strongly on the 
damping of the dynamics of if e xt by the impedance that 
is shunting the SQUID. Experimental control over the 
damping, requires the fabrication of a shunt circuit in the 
direct vicinity of the SQUID, such that its impedance is 
well defined up to the frequency of the SQUID's plasma 
oscillations (microwave frequencies). The shunt circuit is 
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therefore preferably realized on-chip (Z s h in Fig. [I]a). The 
escape from the well may be thermally activated, but 
for underdamped systems with low-capacitance junctions 
quantum tunneling through the barrier can dominate the 
escape rate at low temperatures. The influence of the 
damping circuitry on the Isw statistics f| [f3p7| ] is now 
well understood. A SQUID with very underdamped dy- 
namics usually has Isw values much below Ic, and his- 
tograms of a set of Isw recordings will be very wide. This 
means that one needs to average over many repeated mea- 
surements to achieve the required resolution in readout. 
Thereby, averaging also needs to take place over many 
repeated experiments on the qubit, such that only a time- 
ensemble average can be measured. With a shunt that pro- 
vides high damping at the plasma frequency very narrow 
switching current histograms can be realized @|||i^,[44|,[45| 
that in principle allow for single-shot readout in qubit ex- 
periments. While in such a scheme the SQUID's noise will 
also be enhanced, it is possible to engineer (for realistic 
fabrication parameters) a shunt impedance that is at the 
same time compatible with coherent dynamics of the qubit 
and single-shot readout [[38|. The engineering of single- 
shot readout will not be addressed in detail in this paper. 

The main disadvantage of the switching SQUID is that 
it is not very efficient. During each cycle through the hys- 
teretic TV it is only measuring for a short time. More- 
over, the TV is very nonlinear, such that the repetition 
frequency must be an order lower than the bandwidth of 
the filters. The filtering that is required for realizing low 
effective temperatures and the SQUID's shunt circuit have 
typically a bandwidth well below 1GHz, and the accurate 
readout electronics set a similar limit to the bandwidth. 
In practice this limits the repetition frequency to values in 
the range of 10 kHz |0,||1 to 1 MHz More efficient 

readout may be realized with AC readout techniques (see 
e. g. Ref. ||). 

The slow operation of the switching DC-SQUID sets 
requirements for the mixing rate T r of the qubit. It needs 
be longer than the time required to perform a switching 
current measurement, which requires a time in the range 1 
/us to 100 fis. One could go to shorter times by setting the 
SQUID ready at a high bias current when an experiments 
on the qubit is started, but it is also needed to have the 
mixing time longer than the time it takes to ramp the 
bias current through the range of the switching current 
histogram. At the same time we should realize that the 
quantum system is prepared by waiting for it to relax to 
the ground state, so relaxation times very much longer 
than 100 fis will prohibit a high repetition frequency. A 
high repetition frequency is needed if the signal can only 
be build up by averaging over many switching events. 

The experiments aim at working with many coherent 
Rabi oscillations with a period of about 10 ns [jl9). We 
therefore aim at engineering SQUIDs that cause a dephas- 
ing time that is much longer than 10 ns. The dephasing 
and relaxation times turn out to be shortest at high bias 
currents through the SQUID. Unless mentioned otherwise, 
we will make in this article worst case estimates for the 



dephasing and relaxation times using bias current values 
near the switching current. 

3.2 On-chip control lines 

An attractive feature of macroscopic qubits is that one 
can address individual qubits with control signals from 
microfabricated lines (see also Fig. |l|b,c). For persistent- 
current qubits, for example, a supercurrent through a line 
that is mainly coupled to one specific qubit can be used 
for tuning this qubit's energy bias e. Also, it is convenient 
to provide the microwave signals for control of the qubit's 
quantum dynamics using local superconducting lines. If 
this is realized in a microwave cavity with its first res- 
onance well above the applied microwave frequency, one 
can apply microwave bursts with fast switch times with- 
out being hindered by high-Q electromagnetic modes in 
the volume that is formed by the cold metallic shielding 
that surrounds the sample. 

Microwave signals can be applied using external mi- 
crowave sources at room temperature. Alternatively, on- 
chip oscillators for example based on Josephson junction 
circuits J48| , [l9| can be applied. High microwave currents in 
the control lines are achieved by shorting the microwave 
coax or wave guide close to the qubit with an inductance 
that has an impedance much lower than the source's out- 
put impedance (Fig. [l]b) . For external microwave sources, 
the typical level for the output impedance will be that 
of the available coax technology, typically 50 Q. With on 
chip Josephson oscillators the typical output impedance 
is one order lower. In both cases, it is in practice very te- 
dious to engineer these impedance levels and our analysis 
below will show that this forms a constraint for qubit ex- 
periments: long decoherence times are in conflict with the 
wish for local qubit control and low power levels of the 
applied microwave signals. 

If one uses external microwave sources at room temper- 
ature it is harder than for the quasi DC signals to filter 
out the high temperature noise. Low effective tempera- 
tures can be achieved by a combination of narrow-band 
microwave filters and strong attenuators at low tempera- 
tures. 

4 Relaxation and dephasing from a switching 
DC-SQUID 

4.1 Current-phase relations for the DC-SQUID 

The DC-SQUID has two phase degrees of freedom, the 
gauge-invariant phases 7 r and 7/ of the junctions f}6|| . 
They are related to the supercurrents through the left and 
the right junction, 

h = (Ico + # a )sin7 i , , . 

I r = (Ico-^)sm lr . (0> 

Here I co is the average of the critical current of the two 
junctions. A small asymmetry in the junctions' critical 
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currents is accounted for by AI co <C I co (typically a few 
percent). We will work here with the sum and difference 
phase coordinates Pint and p ex t, which are related by a 
linear transformation 



Pext 
Pint 



7i+7r 

2 



72 = Pext + Pint 
7r = Pext — Pint 



(16) 



The new phase coordinates are related with the current 
passing through the SQUID I sq and the circulating current 
in the SQUID I cir 



Isq — h + Ir 
± ctr — 9 



i =-i —I 

± r 2 S( i ± cir 



(17) 



yielding the following current-phase relation for I sq and 

Icir 



I sg = 2I co COS (pint sin ip e 



AI co sin (p int cos (p ex ^18) 



Icir — Ico 

sin (p^t cos ip ex t + -zAI co cos (p int sin p ex (19) 

We will assume that the DC-SQUID has junctions with a 
critical current and capacitance that are lower than that 
of the qubit junctions. In this case, the internal phase (pint 
follows the flux adiabatically up to time scales much faster 
than -. We will therefore use 



<2> del , 
Pint = 7T— = /• 



(20) 



4.2 Noise on the qubit from the DC-SQUID resulting 

in J((jj) 

The noise that is induced by the measuring SQUID results 
from Johnson- Nyquist noise of the total impedance Z t (u>) 
between the leads that are attached to the SQUID. The 
impedance Z t (ui) is formed the SQUID 's impedance in 
parallel with the impedance of the wiring and circuitry 
that the SQUID is connected to (see the circuit models 
in Fig. |J). At bias currents well below the critical current 
Ic, the phase dynamics can be linearized and the SQUID 
can be modeled as an inductor Lj. The coupling of (p ex t 
to the SQUID's inner degree of freedom ipi n t and thereby 
to the qubit slightly alter the effective value for Lj, but 
the correction it is so small that it can be neglected. The 
Fourier-transformed power spectrum (SV(t) <5U(0)) W of 
the Johnson- Nyquist voltage fluctuations SV across the 
SQUID is MM 



(SV SV) U = huRe{Z t (uj)} coth 



2k B T 



(21) 



We will now calculate how this voltage noise leads to fluc- 
tuations Se of the energy bias on the qubit. As a rule, 
the spectral density J(co) in (g) can then be derived by 

dividing the expression for (Se Se)^ by ti 2 coth ^fc^r ) ■ 

The current-phase relations for I sq and I C ir can be used 
for expressing the current fluctuations. The first term of 
@ gives 



disc, ■ T 



2I co cos / cos (p 

2I co COS /COS Pcxt^V, 



ext dt 



(22) 



where we used (p e xt for the time average of <p ex t- With a 
similar expression for the second term of ( |l8| ) the current 
fluctuations in I sq are 

5I sq w (21 co cos / cos (pext - AI co sin / sin (p ext ) 5(p ex t- 

(23) 

The SQUID is usually operated in regions where the av- 
erage external flux in its loop is between an integer and 
half-integer number of ^o- At these points |cos /| |sin / |. 
Therefore, the second term in ( p3| ) can be neglected unless 
\I co cos (pext\ ~ | AI co sin (p ext | . That is, it can be neglected 
unless the bias current is very high, for which sin Cp ex t ap- 
proaches 1. For most purposes we can thus use 



3Isq ~ 2I co COS / COS (p e xt Sep 

ext ■ 

This is also used to define Lj by expressing 

y _ , dlsn 



dt 



(24) 



(25) 



such that with (g2j), (|20|) and (|18|) Lj should be defined 
as 



Lj = 



1 



1 



2e 2I co cos / cos p e 



2c 



4J2 cos 2 /-4 2 , 



sq 



(26) 



For I c ir we get a similar expression as ( f23[ ) 

Sldr ~ (-/ co sin/sin(^ e2;t + ^AI C0 cos / cos (p ext ) Sp ext . 

(27) 

Using again that the SQUID is operated where |cos/| « 
| sin /| shows that the second term in ( fj~7| ) can be neglected 
unless \I co siri(p ex t \ ^ | AI co cos (p ex t \ ■ For SI C i r the second 
term only plays a role at low bias currents in the SQUID 
for which (p e xt ~ 0, and for most purposes we can use 



Sir. 



Ico 

sin / sin ip ext Sp 

e 



(28) 



In the above we used (f ex t for the time average of <p ex t, 
but at places where it is not confusing it will be simply 
denoted as (p e xt- 

Both noise in I sq and I C i r can couple to the qubit, but 
we will assume that the qubit is mainly sensitive to noise 
in I c ir (as in the experiments in [T^] , where the qubit was 
placed symmetrically inside the SQUID's loop) and ne- 
glect an inductive coupling to noise in I sq . For a more 
general approach, coupling to noise in I sq can be treated 
on a similar footing as noise in Idr-, but for all useful sam- 
ple geometries it should give a contribution to relaxation 
and dephasing rates that is at most on the same order as 
that of Idr- 

With iuj8I C ir — — ^ico sin/ sin </3 e:Et SV follows for the 
fluctuations 5I C i r 



(SIcir SI C 



^) \l 2 co sin 2 / sin 2 pext (SV SV) L 



(29) 



The fluctuations in the imposed qubit flux are S<P q = 
MS I C ir, where M the mutual inductance between the SQUID 
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loop and the qubit loop. This then yields the fluctuations 
in the energy bias with Se = 2I p 5<P q: 



(Se Se) L 



4 



2 M 2 I 2 p I 2 co sin 2 / sin 2 Vext (SV SV) U 



where I v is the amplitude of the circulating current in the 
qubit in the semiclassical states. Using ( |l8| ) and ( |2l| ) and 
filling in = @o this can be written as 



(Se fe)^ft(2^) 2 l(^) 2 / 2 9 tan 2 / 
3 {Z t (a;)}coth( 5 ^) 



Rc| 



(31) 



The fluctuations (Se Se) w are the result of the coupling 
to the oscillator bath, as in (Q). This can be used to define 
J(oj) for our specific environment, 

JH = (^) 2 Is q ^fRe{Z t (u)}. (32) 

These results show, that although the SQUID is perma- 
nently close to the qubit, it may be effectively decoupled if 
there is no net bias current I sq flowing through the device. 
The physical reason for this becomes apparent in Eqs. ( p"8| ) 
and (|l9|): The SQUID remains mirror symmetric in that 
case and consequently the fluctuations of the bias current 
are diverted symmetrically around the arms of the SQUID 
and do not produce flux noise f5l[] . 



4.3 Relaxation times 

With (|) and (||) follows the SQUID 's contribution to 
the relaxation rate. It is here expressed as a function of 
the resonance frequency uj res = v/h at which the qubit is 
operated, 



r r = 



(32) 



2ft 

Re{Z t (uj rt 



/ s 2 9 tan' 



fx 



,)}coth( 



2k B T 



(33) 



In this formula one can recognize a dimcnsionlcss factor 
(~l^f ) which is a scale for how strongly the qubit is cou- 
pled to the measuring SQUID. A dissipation factor in the 
form I 2 R can be recognized in J 2 9 tan 2 / Re{Z t (uj)}. The 
dissipation scales with the absolute value of the current 
fluctuations, so with I sq , and the expression is indepen- 
dent of the critical current of the SQUID junctions I co 
(unless Re{Z t (uj)} depends on I co ). A weak measurement 
scheme in which the inductive coupling to a DC-SQUID 
(M/p/^o) 2 << 1 can yield relaxation rates that are very 
low when compared to a scheme in which leads are di- 
rectly attached to the loop |52] . A measurement of such a 
scheme's switching current could also be used for probing 
the qubit, but the influence of the voltage noise would be 
dramatically worse. 

With the result (^3|) the relaxation rate for typical 
sample parameters will be calculated. Sample parameters 



similar to our recent experiment 1141 are io res = 10 GHz, 
A = 2 GHz, ^ = 0.002. It is assumed that a SQUID 
with 21 co = 200 nA is operated at / = 0.75 ir and bi- 
ased near the switching current, at I sq — 120 nA. For 
T = 30 mK the relaxation rate per Ohm environmental 
impedance is then 



r: 



150 fisn 

Re{Z t (oj res )}' 



(34) 



4.4 Engineering Rc{Z t {uj)} for slow relaxation 

In practice the SQUID's resolution is improved by building 
an on-chip electromagnetic environment. We will consider 
here a large superconducting capacitive shunt (Fig. |]a, as 
in our recent experiment |14| ). This scheme will be de- 
noted as the C-shunted SQUID. As an alternative we will 
consider a shunt that is a series combination of a large ca- 
pacitor and a resistor (Fig. [|b). This will be denoted as the 
i?C-shunted SQUID. The C shunt only makes the effec- 
tive mass of the SQUID's external phase tp ex t very heavy. 
The RC shunt also adds damping at the plasma frequency 
of the SQUID, which is needed for realizing a high reso- 
lution of the SQUID readout (i. e. for narrow switching- 
current histograms) The total impedance Z t (uj) of 
the two measurement circuits are modeled as in Fig. ||. 
We assume a perfect current source Ibias that ramps the 
current through the SQUID. The fact that the current 
source is non-ideal, and that the wiring to the SQUID 
chip has an impedance is all modeled by the impedance 
Zi. The wiring can be engineered such that for a very 
wide frequency range the impedance Z\ is on the order of 
the vacuum impedance, and can be modeled by its real 
part R[. It typically has a value of 100 fl. On chip, the 
impedance is formed by the Josephson inductance Lj in 
parallel with the shunt circuit (C s h, or the series combi- 
nation of R s h and C s h)- We thus assume that the total 
impedance Z t (uj) can be described as 




(35) 



where R s h should be taken zero for the C-shunt scenario. 
The circuits in Fi g. H are damped LC resonators. It 
(0) and Iffi ) that one should keep the LC- 

where Re{Z t (w)} 



is clear from ([7|) and 
resonance frequency u>lc 



has a maximum, away from the qubit 's resonance u) r es = 
v/h. For practical values this requires lo^c <C Lo res for Alu- 
minum technology (with Niobium-based technology, the 
regime ljlc ^ w res may be realized [0). This then gives 
the circuits a Re{Z t (u;)} and J(ui) as plotted in Fig. [|. 
For the circuit with the C shunt 



Re{Z t (w)} 



, , 2 T 2 

Ri ' 

Ri, 



for oj <C ojlc 
for lo = U)lc 
, for iv > lo lc 



(36) 
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For the circuit with the RC shunt 

uj 2 L~ 



term in (||) dominates, such that with ( f4l| ) and for the 
qubit operated where e as v 



Re{Z t (uj)} 



Ri , for lo < lolc 
<Rl, for u = lolc < r^c^ 



Rl//R s h, for w 
Ri//R s h, for w » Wic 



(37) 



The difference mainly concerns frequencies lo > lolc, where 
the C-shunted circuit has a Re{Z t (uj)}, and thereby a re- 
laxation rate, that is several orders lower than for the 
-RC-shunted circuit. 

For a C-shunted circuit with lolc <C Lo res the Re{Zt(LO re s 
— 2 — 772— 5-- This yields for J(w) at lo > lolc 



J{uo) 



(2tt) 1 /M/j 



#0 



7 S 2 tan 2 / 



(38) 



The factor indicates a natural cut-off for J(lo), which 
prevents the ultraviolet divergence J|(],[32| and which in 
much of the theoretical literature is introduced by hand. 
The i?C-shunted circuit has softer cut off 1/lo. The mixing 
rate for the C-shunted circuit is then 



(A/%f (2tt)- 



2h 



1% tan 2 



/ 



coth 



( hu ra , \ 

\2k B f J ■ 



(39) 



Fig. g presents mixing times r r vs Lo res for typical sample 
parameters (here calculated with the non-approximated 
version of Re{Z t (uj)}). With the C-shunted circuit it seems 
possible to get r r values that are very long. They are com- 
patible with the ramp times of the SQUID, but too slow 
for fast repetition rates. In Fig. one can directly see 
from the values of J(lo) that an i?C-shunted circuit with 
otherwise similar parameters yields at Lo res /2n = 10 GHz 
relaxation times that are about four orders shorter. For 
the parameters used here they are in the range of 15 /is. 
While this value is close to the desired order of magnitude, 
one has to be aware of the fact at these high switching cur- 
rent values the linearization Eq. (26) may underestimate 
the actual noise. In that regime, phase diffusion between 
different minima of the washboard potential also becomes 
relevant and changes the noise properties |37 53 . 



4.5 Dephasing times 

At low frequencies lo < lolc the C-shunted and RC- 
shunted scheme have Re{Z t (u>)} w R J such that with 
© and (|) 



J(«) 



(27T) 2 



2tt 

T 



MI, 
<P 



M ^tan 2 /^, 



MIp 



(40) 



(41) 



The environment is Ohmic at low frequencies since we 
have J(lo) cx lo. For our sample parameters the second 



(2tt) 2 ( MIp 
h 2 \ <P 



I sq tan 



2 / f k B T. 
Rl 



(42) 



Note in ( p6[ ) that Lj oc l/I co , such that the dephasing rate 
( |42"| ) does not depend on the absolute value of the current, 
but on the ratio I sq /Ic- For the typical sample parameters 
as used in Fig. || the dephasing time is about 10 ns, which 
is too short. However, we can gain a few orders (if r r is low 
^ngugh) by the fact that we can do the quantum coherent 
control at low I sq (the previous estimate was calculated 
for I sq — 120 nA, in the switching region). At I sq — we 
find r$ = in this linear approximation for the SQUID 
inductance. At I sq = we should therefore estimate the 
dephasing due to second order terms. However, in practice 
the dephasing is probably dominated by the second term 
in (p7j), which is due to a small asymmetry in the fabri- 
cated SQUID junctions of a few percent. This influence 
can be mapped on a small bias current (a few percent of 
the critical current, say 5 nA) through the SQUID. There- 
fore, at I sq w the dephasing times can be ("Hp) 2 times 
longer. Furthermore, the factor L 2 , as defined in (|26|), is 
at 5 nA about a factor 2 lower than at 120 nA. For our 
parameters this allows for w 20 //s, see also Fig. |?]. Fur- 
ther improvement is feasible by making e. g. Ri = 1 ki7, 
working with a lower mutual inductance M or tuning the 
qubit close to the degeneracy point, as in |]l2| . 

Finally, we would like to mention that in the literature 
on dissipative two-level systems one often assumes Ohmic 
dissipation, corresponding to a purely resistive shunt across 
the junctions of the qubit. For a description of such a sys- 
tem one usually introduces an artificial exponential cut off 
at frequency lo c , yielding J(lo) of the form 



J{lo) = acjcxp( ). 

LOr 



(43) 



In our case, J{lo) has substantial internal structure orig- 
inating from the frequency-dependence of Re{Z t (Lo)} . In 
order to compare our results to the Ohmic case, we plot in 
Fig. | an Ohmic fit to the actual J(co) of the DC SQUID. 
For the parameters as in Fig. the resemblance is reason- 
able for a resistive shunt corresponding to a = 0.00062, 
and a cut off lo c = 0.5 GHz. For low currents, as for the 
dashed line in Fig. a = 1 • 10~ 7 . This corresponds to an 
extremely underdamped system, with a long dephasing 
time. 



5 Relaxation and dephasing from on-chip 
control circuits 

We will treat here the influence of noise from the mi- 
crowave leads in a similar way as worked out for the SQUID. 
Here the total environmental impedance Z t (Lo) is formed 
by a 50 fi coax, that is shorted at the end by a small 
inductance L mw , see the circuit model in Fig. pi This in- 
ductance L mw has a mutual inductance M mw to the qubit. 



10 



Caspar H. van der Wal et at: Engineering decoherence in Josephson persistent-current qubits 



The voltage noise is given by (|T]) . The noise leads to fluc- 
tuations 8s of the energy bias separation s as follows. The 
current noise in L mw is 81 l — j^jj; — 8V. The qubit fluc- 
tuations in the flux <l> q are then 8<P q = M mw 8lL, and the 
fluctuations in the energy bias are 8e = 2I p 8'P q . This gives 
for the fluctuations (8s Ss)^ 



(8s 8s) u 



4ft / M mw I p 



Re{Z t (w)}coth 



2k B T 



and for the environmental spectral density (||) 



(44) 



JH= 4 fM^lA Re{ZtH} (45) 

lUv \ ^mw J 

The Re{Z t (ui)} is that of a first order low-pass LR filter 
with a -3 dB frequency w^jj = R/L. For the L mw to be ef- 
fectively a short its impedance u>L mw should be small com- 
pared to 50 fl at the frequency of the applied microwave 
radiation (typically 10 GHz), giving L mw <C InH. This 
can be realized by making the length of the short line less 
than about 100 /xm. This means that all relevant frequen- 
cies are below the -3 dB frequency ll>lr, and that for both 
relaxation and dephasing we can approximate 



Re{Z t H}«^^ 



(46) 



with R mw = 50 fl. For oj < lolr, Ji 1 ^) is again Ohmic, 



J{w) 



4u (M mw I P Y 



Rn 



and with 



we find for a 



4 (M mw I P Y 
2irh Rmw 



(47) 



(48) 



Note that these results are independent of L mw . A larger 
L mw leads to enhanced voltage noise, but the resulting 
current noise is reduced by the same factor. For frequen- 
cies below ujlr the current noise is just that of a shorted 50 
fl resistor. For frequencies higher than lolr, Be{Z t (w)} w 
Rmw, such that J(uj) has a very soft intrinsic l/u> cut off. 

For the relaxation rate (0) as a function of the qubit 's 
resonance frequency uj res we now have 



2A 2 



(M mw I p ) 



Rn 



■ coth 



\2k B T 



(49) 



This has a much weaker dependence on uj res than for the 
SQUID, results are plotted in Fig. || The results are plot- 
ted for M mw — 0.1 pH (further parameters are as used 
for the SQUID calculations) and for this M mw value the 
relaxation times are in the required range of about 100 
[is. The value M mw 0.1 pH corresponds to a 5 (im loop 
at about 25 /im distance from the microwave line, and is 
compatible with the fabrication possibilities and the mi- 
crowave requirements. With this geometry it is still possi- 
ble to apply sufficient microwave power for pumping the 



qubit 's Rabi dynamics at 100 MHz 
an oscillating <P q of about 0.001 



Lz (L 
■o @ 



e. pumping with 
, which needs an 



oscillating current of ss 20 /iA, corresponding to 

20 nW, i. e. — 47dBm microwave power), while the dis- 
sipated microwave power in the attenuators at the refrig- 
erator's base temperature remains well below the typical 
cooling power of 1 fiW. 

For the second term of the dephasing rate (B) we thus 
find for the qubit operated where s ~ v 



^ 4 (M mw I p ) 
ft? R B 



(50) 



Using the same parameters as in the above calculation of 
the relaxation time we find awl- 10~ 7 and for T = 30 mK 
the dephasing time is w 130 /is. While this dephas- 
ing rate is sufficient for demonstration experiments and 
promising for applications, we like to note that it is much 
harder to engineer this dephasing rate as compared to the 
DC-SQUID dephasing. It is in practice quite tedious to 
apply microwave technology with impedance levels R m w 
much higher than 50 S7, both for externally generated 
microwaves and on chip generators (it could for instance 
be increased using a planar impedance transformer flSjJ). 
Making M mw smaller requires higher microwave currents, 
and thereby more microwave dissipation on the mixing 
chamber. The cooling power per qubit will quite likely re- 
main below 1 fJ.W, so much stronger microwave signals 
from a larger distance is not an option. Moreover, mak- 
ing M mw very small means that the control line is 100/xm 
or further away from the qubit. In this case, applying mi- 
crowaves locally to one specific qubit on a chip with several 
coupled qubits is much harder. 



6 Suppressing rates by freezing states and 
idle states 

With additional control techniques the decoherence rates 
can be made better than the estimates made in the pre- 
vious sections. These are based on the pre-factors in Eqs. 
(0) and (H) Bringing the qubit in a so-called idle state 
(e = 0) [|L8| can reduce the dephasing, but not beyond 
r r /2 (which is enhanced at e = 0). 

Another useful technique is freezing. Here, the tun- 
nel coupling A between the two qubit states is strongly 
reduced before the measurement process starts with a 
fast but adiabatic control current [^9). This allows for 
much slower measurements, thus for weaker coupling to 
a damped SQUID with very high resolution. Moreover, 
it has the advantage that the tunnel coupling becomes 
so week that the Hamiltonian almost commutes with a z . 
This can improve the correlation between the outcome of 
a z measurements, and energy states in the case that the 
calculation states are energy eigen states. However, freez- 
ing requires that the qubit junctions are realized as small 
DC-SQUIDs. This means that a x noise will be strongly en- 
hanced. The influence of this additional noise source can 
be calculated along the same lines as in (32) and, assum- 
ing that it is not correlated with the other noise sources, 
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add up to the rates we have calculated in this paper. See 
Fig. H for numerical estimates on the possible adiabatic 
suppression of A 



7 Discussion and conclusion 

We have developed a scheme for modeling the decoherence 
of a persistent-current qubit due to its electromagnetic en- 
vironment. Examples for both control and read-out elec- 
tronics are worked out quantitatively. We discussed how 
the dephasing and relaxation rates of the qubit can be de- 
rived from the impedance of the electronic circuitry, and 
provided design criteria for such electronics. In particular, 
we have shown that even though the readout SQUID is 
always close to the qubit, it can be effectively decoupled. 
Our examples show that the present status of experimen- 
tal technology should allow for the observation of quan- 
tum coherent oscillations between macroscopic persistent- 
current states. 

In this final section, we like to point out that the theory 
that is used is still in development. In particular, an envi- 
ronment with a strongly structured spectrum may violate 
the weak-coupling Born approximation at its resonance 
frequencies, but it may also induce weak additional deco- 
herence off-resonance. This situation is hardly accessible 
with traditional theoretical methods for this problem and 
alternative approaches such as the flow-equation scheme 
in [p5| may be needed. Furthermore, it is not clear whether 
the impedance can be described by a single temperature. 
At low frequencies, noise from parts in the system with 
a higher temperature can reach the sample. However, the 
experimental results reported in Ref. |12| indicate, that 
an analysis of the decoherence of the type we give here 
gives good predictions of the experimental decoherence 
time scales in a superconducting qubit. 
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Fig. 1. Experimental setup for measurements on a Josephson 
persistent-current qubit. The qubit (center) is a superconduct- 
ing loop that contains three Josephson junctions. It is induc- 
tively coupled to a DC-SQUID (a), and superconducting con- 
trol lines for applying magnetic fields at microwave frequencies 
(b) and static magnetic fields (c). The DC-SQUID is realized 
with an on-chip shunt circuit with impedance Z s h- The circuits 
a)-c) are connected to filtering and electronics (not drawn). 
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Fig. 2. Sketch of a typical hysteretic current- voltage 
characteristic (IV) for a current-biased Josephson junc- 
tion or small DC-SQUID. The IV is hysteretic; arrows 
indicate which of the two branches is followed at an in- 
crease or decrease of the bias current. When the bias 
current Iuas is ramped up from zero (a), the voltage V 
first remains zero. The circuit is here on the supercurrent 
branch of the IV (b). When Ibias approaches the criti- 
cal current Ic, a slow diffusive motion of the phase tp e xt 
leads to a very small voltage across the system (c). At 
slightly higher current (d), but always below Ic (e), the 
system switches to a running mode for (p ex t, and the volt- 
age jumps to a value set by quasiparticle tunneling over 
the superconducting gap, V = 2A/e (this current level 
(d) is the switching current Isw)- At further increase 
of the current (f) the IV approaches an Ohmic branch, 
where transport is dominated by quasiparticle tunneling 
through the normal tunnel resistance of the junctions. 
When lowering the bias current the system follows the 
running mode (g) down to a low bias current where it 
retraps on the supercurrent branch (at the level I re trap, 
indicated by (h)). See also the corresponding washboard 
potential model, in Fig. H. 




Fig. 3. The dynamics of a current-biased DC-SQUID, 
modeled as a particle with coordinate <p ex t in a one- 
dimensional tilted washboard potential U. The labeling 
a)-f ) corresponds to that of Fig. ^. At zero (a) and small 
bias currents (b), the particle is trapped in a well of the 
washboard. Apart from the small plasma oscillations at 
the bottom of the well, the particle's coordinate <p ex t 
is fixed. When increasing the slope of the washboard, 
the particle will start to have a slow, on average down- 
wards, diffusive dynamics, with rare excursions to one of 
the neighboring wells (c). At the switching current Isw 
there is a high probability that the trapped particle will 
escape to a running mode (d) , with effectively zero prob- 
ability for retrapping. Here the loss of potential energy 
exceeds the dissipation when the particle moves one pe- 
riod down the washboard, and the particle builds up a 
high kinetic energy. Due to thermal fluctuations, exter- 
nal noise, and in certain cases quantum fluctuations, this 
occurs below the critical current Ic'- the slope where all 
local minimums in the washboard potential disappear 
(e). At currents higher than this slope (f), the particle 
will always be in a running mode. The retrapping process 
when lowering the bias current follows similar dynamics. 
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Fig. 4. Circuit models for the C-shunted DC-SQUID 
(a) and the i?C-shunted DC-SQUID (b). The SQUID 
is modeled as an inductance Lj. A shunt circuit, the 
superconducting capacitor C s ;, or the R s h-C a h series, is 
fabricated on chip very close to the SQUID. The noise 
that couples to the qubit results from Johnson-Nyquist 
voltage noise SV from the circuit's total impedance Zt. 
Zt is formed by a parallel combination of the impedances 
of the leads Zi, the shunt and the SQUID, such that Z t = 
(1/Zi + l/(R sh + l/iuC sh ) + 1/iuLj-y 1 , with R sh = 
for (a). 



C shunt 



RC shunt 
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Fig. 5. A typical Re{Z t (uj)} for the C-shunted SQUID 
(a) and the _RC-shunted SQUID (b), and correspond- 
ing J(u>) in (c) and (d) respectively. For comparison, 
the dashed line in (c) shows a simple Ohmic spectrum 
( p3| ) with exponential cut off lo c /2ty — 0.5 GHz and a 
= 0.00062. The parameters used here are I p = 500 nA 
and T = 30 mK. The SQUID with 21 co = 200 nA is 
operated at / = 0.75 7r and current biased at 120 nA, a 
typical value for switching of the C-shunted circuit (the 
i?C-shunted circuit switches at higher current values). 
The mutual inductance M = 8 pH (i. e. ^ = 0.002). 
The shunt is C s h = 30 pF and for the RC shunt R s h = 
10 n. The leads are modeled by R t = 100 Q. 




co 12% (GHz) 



Fig. 6. Typical relaxation times due to the C-shunted 
SQUID, the 7?C-shunted SQUID, and coupling the mi- 
crowave leads as a function of the resonance frequency 
at which the qubit is operated. The example of the mi- 
crowave leads contribution is for a mutual inductance 
Mmw to the coaxial line of M mw =0.1 pH. Parameters 
are further as described in the caption of Fig. || 




Fig. 7. The dephasing time (^2|) as a function of the bias 
current I aq through the SQUID (solid line). The dashed 
line shows for I aq — 5 nA, a typical minimum value 
for the effective bias current for a SQUID with a few 
percent asymmetry between its junctions. At this point 
i~0 = 131 fj,s, and a = 1 ■ 10~ 7 . Parameters are further as 
described in the caption of Fig. |^. 
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Fig. 8. Circuit model for coaxial line that is inductively 
coupled to the qubit. The coaxial line is modeled as a 
50 O impedance that is shorted near the qubit with an 
inductance L mw . The qubit is coupled to this short with 
a mutual inductance M mw . The noise that couples to 
the qubit results from Johnson-Nyquist voltage noise SV 
from the circuit's total impedance Zt, formed by a par- 
allel combination of the 50 Q impedance and L mw . 
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Fig. 9. Numerically simulated suppression of A in qubits 
where one (solid) or three (dashed) of the junctions are 
realized as small DC-SQUIDs, using typical qubit pa- 
rameters as mentioned in the text. The horizontal axis 
is the factor by which the Josephson energy of the rel- 
evant junctions is increased. For the solid line, one can 
observe a non-monotonic dependence of A on the freeze 
factor, which results from the fact that if the Joseph- 
son coupling of the weak junction is strongly increased, 
alternative tunneling paths between the current states 
open up, see also Ref. Il9|. 



